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Abstract 

An exact solution of tlie Maxwell equations in Rindler coordinates 
is obtained. The electromagnetic field represents a wave preserving its 
shape in a relativistic uniformly accelerated frame. The relation with 
Airy beams is shown explicitly in the non-relativistic limit. 

PACS: 41.20.Jb;42.25.Bs 
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1 Introduction 



In 1979, Berry and Balazs [Ij obtained a solution of the Schrodinger equation, 
in terms of an Airy function, describing a non-spreading accelerating wave 
packet. More recently, the Airy beam, an equivalent solution of the paraxial 
equation of optics, has received a considerable amount of attention both from 
a theoretical as well as from an experimental point of view (see, e.g., Ref. 

These previous works on accelerating waves have been restricted to the 
non-relativistic case. The aim of the present work is to obtain an exact 
solution of the Maxwell equations in Rindler coordinates, which are the nat- 
ural coordinates associated with uniform acceleration in special relativity. 
The ensuing exact solution represents an electromagnetic wave that is shape 
invariant in a uniformly accelerated frame and appears therefore as an accel- 
erating wave in an inertial frame. The relation with the Airy beam is shown 
explicitly in the non-relativistic limit. 

In Section 2, the Maxwell equations are solved in Rindler coordinates 
for the two fundamental modes of the field. A particular solution, invariant 
under space-time hyperbolic rotations, is obtained and the explicit forms of 
the energy density and the Poynting vector are worked out. The relation 
with the Airy beam in the non-relativistic limit is shown in Section 3 and a 
brief discussion of the results is given in Section 4. 



2 Maxwell equations 

The Rinder metric is obtained with a change from Cartesian {ct, x, y, z) to 
Rindler coordinates (cT, x, y, Z) defined as 

ct = Z sinh (aT), 

z = Zcosh(aT), (2.1) 

where a is a parameter with dimensions of frequency (or acceleration divided 
by c). The metric thus takes the form 

ds^ = -a^Z^dT^ + dx^ + dy^ + dZ\ (2.2) 

The electromagnetic field tensor Fap follows from a vector potential Aa 
such that 

Fai3 = do,Ap - dpA^ (2.3) 
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and is assumed to satisfy the source-free Maxwell equations 

d 



fzF'^n = 0. (2.4) 



The two modes of the electromagnetic field in Rindler coordinates can be 
represented by the transverse electric (E) and transverse magnetic (M) vector 
potentials, 

Al = (0; a2t/^, - aiC/^, 0) (2.5) 

and 

>lf = (aZSaC/^; 0, 0, ^^oC/^) , (2.6) 

given in terms of two potentials and . The total electromagnetic field 
turns out to be 

Fi2 = -dl^U^ - dljj"^ , (2.7) 
aZ 

Foi =-aZdl^U'^ + dl^U'' , 
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/ aZ 



As can be checked by direct substitution, the complete Maxwell equations 
are satisfied if 

d^oU - a^Z^ VlU - c?Z {z d^U^ = , (2.8) 

for both potential and . Here Vi = dl^ + dl^. 

This last equation admits a solution by separation of variables: 

U{xi') = e-'^'^+'''-^+'''yyu{Z), (2.9) 

where Q is the frequency as measured in the accelerated frame and U{Z) 
satisfies the equation 



(-V-klZ'h^O, (2.10) 

\ a ^ 



for both modes, with fc^ = k1 + kf,. The solution is 

U{Z) = C^E,M)K,n/a{k±Z), 



(2.11) 



where Ki^/a{k±Z) is the modified Bessel function of the third kind with 
imaginary order and C(e,m) are constants. The other hnearly independent 
solution is Lifi/a{k±Z), but it diverges at Z — >■ oo [5]. 

The electromagnetic tensor has the following forms in Rindler coordi- 
nates. For the transverse electric mode: 



„—iQT+ikxX+ikyy 



{azy 







c-^nkyK 



-c'^nkyK 







c~^Qk^K 
kliaZfK 
c-^VLk^K -kliaZyX 

ikyk^{aZfK' -ik^k^{aZfK' 



\ 

-ikyk^{aZfK' 
ik^k^{aZfK' 

/ 
(2.12) 



and for the transverse magnetic mode: 



M 



-iQT+ikxX+ikyy 



aZ 



( 

ikxkxK' 
ikyk I K' 



-ik^k±K' -ikyk±K' -k^K \ 

-c~^VLk^K 
-c-^VLkyK 

J 













(2.13) 



\ klK c-^Vtk^K c-^VtkyK 

where K = KiQ/a{k±Z) and K' is the derivative of K with respect to its 
argument. 



2.1 Minkowski coordinates 

Transforming to Minkowski coordinates (the electric and magnetic fields are 
Ea = Fao and Ba = \ eabcF''^), it follows that 



E = e 



-iUT+ikxX+ikyy 



CeG X + F k I -P,K e 



g ^—iQT+ikxX+ikyy 

where 



F = — smh.(aT)K — ik± cosh.(aT)K' 
aZ 



(2.14) 
(2.15) 

(2.16) 



G = — : cosh.{aT)K — ik± smh.{aT)K . 



(2.17) 



In these formulas, T and Z must be interpreted as functions of t and z, as 
given by Eq. (12. ip . 

It follows from these equations that the invariants of the electromagnetic 
field are 

E ■ B* + c.c. = -{CeCI, + CMC*E)k\_7{k^Z) (2.18) 

and 

E ■ E* - B ■ B* = {\Cm? - \CE?)kl7{k^Z), (2.19) 

where 



F{k^Z) = {K'Y + 



1 - 



.ak±Z 

As for the energy density W and Poynting vector S, it follows that 



(2.20) 



W = E-E* + B-B* 

and 



{\CEr + \CMnki cosh{2aT) g{k±Z)+K' (2.21) 



S = E X B* + c.c. = kl cosh(aT) x 
{{\CE\' + \CM\y[smh{aT) g{k^Z) e, + 2(-^)irV 

2t{CMC*E - CeCI^KK' x e,}, 



where 



(2.22) 
(2.23) 



-ak^Z , 

and = k^/A;_L is the unit vector in the direction of k_L. 

The total energy integrated over the whole Rindler wedge, z > 0, is 
unbounded due to the divergence of the modified Bessel function at the 
horizon z = ±ct. This is related to the nonphysical behavior of the field 
in that limit, since it would correspond to a section of the electromagnetic 
field with infinite acceleration. Nevertheless, it is possible to integrate from 
a certain value of Z, say e, to Z — )■ cxd and obtain a finite value of the energy. 
According to Eq (10. 7p in the Appendix, 

nk^Q 



w dz = {\CEr + \c. 



M\ 



ea sinh(7rfi/a) 



(2.24) 



the integration being taken over the three-dimensional space t = = T and 
z>0. 
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3 Non-relativistic limit 



The above solutions are exact. The non-relativistic limit can be obtained 
defining first an acceleration parameter g = ac and then performing a coor- 
dinate shift 

z ^ z -\ , 

9 

from where it follows that 

+ (3.1) 
9 2 

and 

sinh(aT) ^ h , cosh(aT) = 1 + 0[{gt/c)\ (3.2) 

and of course T ^ t, thus identifying Q as the frequency in an inertial frame 
in which the velocity of the wave is non relativistic. As a next step, one can 
use some approximation for the Bessel function of imaginary order (see, e.g, 
[SIEIE])- However, as shown by Balogh [7], such approximations break down 
for values of the argument close to the order (the so-called turning point) 
and another kind of approximation must be used involving Airy functions. 
For our present purpose it is enough to notice that, in the non-relativistic 
approximation of Eq. (13. ip , Eq (I2.10p reduces to 



(7)'-'^(^ + 2|c)]f/ = 0, (3.3) 



where = z — ^at"^. The solution in terms of the Airy function Ai is 



U = Cem At 



(3.4) 



The particular case k±c = Q corresponds to the turning point [7] where the 
solution in terms of Bessel functions reduces to the Airy beam treated in the 
literature. To obtain the non-relativistic limit, it is only necessary to change 
the function K of the previous section to the Airy function given above, but 
this does not seem to offer a significant advantage for practical calculations. 
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4 Discussion 



An exact solution of Maxwell's equations in Rindler coordinates was ob- 
tained. From the expression of the Poynting vector in the inertial frame, 
given by Eq. f l2.22p . it is seen that it describes an accelerating light beam 
along the z axis with a turning point at time t = 0. It also has a component 
perpendicular to the averaged direction of propagation, which is typical of 
structured light beams (for instance, Bessel beams [8j). 

The above solution is given in terms of a modified Bessel function of 
imaginary order and is a linear superposition of the two modes of the elec- 
tromagnetic field. In the non-relativistic limit, it reduces to an Airy wave 
for a particular value of the frequency. It should be noticed, however, that 
the total energy of the relativistic accelerating field is not bounded due to its 
divergence at the Rindler horizon z = ±ct; this suggests that the Airy wave, 
if properly extended to the exact relativistic limit, will exhibit a similar di- 
vergence. This fact is not unexpected since the energy of an electromagnetic 
field integrated over the whole space is usually divergent in many idealized 
situations, including plane waves, and it is necessary to constrain the field 
to a finite region of space. Moreover, as shown by Boulware 0, singularities 
at the horizon also occur in the related problem of the electromagnetic field 
produced by a uniformly accelerated charge. 

Appendix 

The modified Bessel function of the third kind is defined in general as 



For purely imaginary order, Kia{x) has the following asymptotic behaviors 
(see, e.g. Dunster [5]): near x = 




(0.1) 



Kia{x) = -(^ 



a sinh(a;7r) 
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and for a; — 7- oo 




(0.3) 
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The energy-density is proportional to the function £{x) = G{x) + Kf^{x), 
where Q is defined by 02.231) . The recurrence relation 



X 



implies 

£{x) = Kia-i{x)Kia+i{x) + Kf^ix). 
Using Nicholson's formula (Watson [TUj. §13.72) 



POO 

K^{x)K^{x) = 2 Kf,^^{2xcosht) cosh(/i ± u)t dt, (0.4) 
Jo 



it follows that 



Since Ki{z) ~ z ^ for 2; <^ 1, we finally obtain 

7Ta 



(0.5) 



poo 

£{x) = 2 / [is:o(2xcosht) + ir2(2x cosht)] cos(2at) dt, 
Jo 

and since Kq + K2 = —2K[, 

1= rS{x)dx= r ^^^^^^^^Ki{2ecosht)dt. (0.6) 

Je Jo COsht 



(0.7) 



e sinh(7rQ;) 
as an approximation valid for e ^ 1. 
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